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Abstract 

^ . Purification is a process in which decoherence is partiaUy reversed by 

using several input systems which have been subject to the same noise. 
The purity of the outputs generally increases with the number of input 
systems, and decreases with the number of required output systems. We 
construct the optimal quantum operations for this task, and discuss their 
Qs , asymptotic behaviour as the number of inputs goes to infinity. The rate at 

O^ ' which output systems may be generated depends crucially on the type of 

purity requirement. If one tests the purity of the outputs systems one at 
a time, the rate is infinite: this fidelity may be made to approach 1, while 
at the same time the number of outputs goes to infinity arbitrarily fast. 
On the other hand, if one also requires the correlations between outputs 
to decrease, the rate is zero: if fidelity with the pure product state is to 
r^ ' go to 1, the number of outputs per input goes to zero. However, if only a 

^H, fidelity close to 1 is required, the optimal purifier achieves a positive rate, 

which we compute. 
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Cl ■ 1 Introduction 

A central problem of quantum information processing is to ensure that devices 
which have been designed to perform certain tasks still work well in the presence 
of decoherence, i.e., under the combined influences of inaccurate specifications, 
interaction with further degrees of freedom, and thermal noise. Decoherence 
typically has the effect of producing mixed states out of pure states, so it is 
natural to ask whether the effects of decoherence can be partially undone, by 
processes turning mixed states into purer ones. As in the classical case this is 
impossible for operations working on single systems. However, if many (say N) 
systems are available, all of which were originally prepared in the same unknown 
pure state ct, and subsequently exposed to the same (known) decohering process 
i?* , then an analysis of the combined state may well allow the reconstruction of 
the original pure state. The quality of this reconstruction will increase with N . 
In fact, it should approach perfection as A^ ^ oo: in this limit one can determine 
the decohered state i?,cr to an arbitrary accuracy by statistical measurements. 
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The question is only, whether the knowledge of the full density matrix Rt:a 
admits the reconstruction of ct, i.e., whether the linear operator R^ is invertible. 
Generically, and for sufficiently small decoherence, this is the case. However, 
the operator R^^ is usually not positive, i.e., it takes some density matrices 
into operators with negative eigenvalues. Therefore, it does not correspond 
to a physically realizable apparatus. But it does describe a computation we 
can perform to reconstruct a from the measured (or estimated) density matrix 
p = R^a. 

How well can this reversal of decoherence be done when the number N of 
inputs is given, and finite? The answer depends critically on the way the purifi- 
cation task is set up, and what "figure of merit" we try to optimize. In general, 
the resulting variational problems may be very hard to solve. However, in the 
specific model situation chosen in this paper, the solution is fairly straightfor- 
ward: we take qubit systems, and assume that decoherence is described by a 
depolarizing channel of the form 

i?,CT = Acr + (l-A)-. (1) 

The purifier will be a device T taking a state of N qubits, and turning out 
some number M of qubits, where M may be cither fixed or itself a random 
quantity. In the latter case T is given mathematically by a family Tm of com- 
pletely positive maps, where Tm takes a density matrix of N qubit systems, 
and produces a positive operator on the M qubit space, which is not necessarily 
normalized to unity: the normalization constant wm — tr(rM(p)) is interpreted 
as the probability of getting exactly M outputs from the input state p. Thus 

Em ^m = 1- ^ 

Our aim is to design T to get outputs as close as possible to the uncorrupted 
input state a, and also as many of them as possible. This is reminiscent of 
cloning problems pi 0]. However, in cloning problems the aim is to get many 
copies of the input state to T, which in our case is the mixed state i?*cr, rather 
than the pure state a. In both cases there is clearly a trade-off between the 
quality of the outputs and their number, which is why there are several different 
ways to state the problem. In the sequel we will briefly describe the variants of 
the purification problem, together with the results, which will be shown later in 
the paper. 

1. Maximal fidelity, failure to produce any output admissible. The best fi- 
delity of outputs is clearly achieved, when the weakest possible demands 
are made on the number of outputs. In this case we do not even insist on 
an output every time the device is run, but only on some non-zero proba- 
bility for getting an output. The best achievable fidelity of these outputs 
goes to 1 as A^ ^ oo, but not substantially faster than with the following 
stronger requirement on output numbers. 

2. M — 1 fixed, number M never increased at expense of output purity. This 
is the approach taken by ||3|. At least one output qubit is required, and 
the figure of merit is based on the fidelity of this one qubit. As it turns out 
the optimal device for this problem can just as well produce more outputs 
of the same optimal fidelity, with a certain rate. However, this rate is not 
part of the optimization criterion. 



3. M fixed, purity measured by one-particle restrictions. For fixed M, N, this 
problem is rather similar to 0. However, with the additional parameter AI 
we can discuss better the trade-ofF between rate and quality of outputs. 
Suppose we fix some dependence of the number of outputs M{N) on the 
number of inputs. Do the states still approach ct as A^ ^ oo? Clearly, if 
M{N) increases slowly, e.g., at the rate given by the optimal device from 
H, this will be the case. What may seem surprising at first, however, is 
that no matter how fast M{N) —^ oo, the state of each output qubit still 
approaches the uncorrupted pure state. In this sense, optimal purification 
works with an infinite rate. 

4. M fixed, purity measured by fidelity with respect to ct*^^. The infinite 
rate depends critically on what we use as the quality criterion for outputs. 
Apart from the fidelity of the restrictions of the output state to single 
qubits used in g we could also look at the fidelity of the outputs with 
respect to the M particle pure state a"^^ , thereby taking into account 
also the correlations between different outputs. For fixed M, the difference 
between these two fidelity measures does not seem so great, because one 
can be estimated in terms of the other. However, the estimates are M- 
dependent (see below), and hence for problems involving a limit M -^ cxd 
the fidelity with respect to the combined state may (and does) turn out 
to be a much tighter criterion. In fact, no process with finite rate M/N 
achieves fidelity— *■ 1, and in this sense even optimal purification works 
with zero rate, in sharp contrast to above. On the other hand, for any 
finite fidelity requirement, there is an output rate for an optimized process, 
which is computed below. 

These results will be stated in precise terms in the following Section ||, together 
with the notation needed for that purpose, and graphs of the optimal fidelities 
and rates. The proofs follow in the subsequent sections. Technically they hinge 
on the decomposition theory of tensor product representations of SU(2), and 
this background is provided in Section H. The reason for representation theory 



to enter in such a crucial way is isolated in Section 3.1, where it is shown that 
the optimal devices can be taken to be SU(2)-covariant (do not single out a 
basis in the qubit space). The two basic purifiers, called the "natural purifier" 
(optimal for question above) , and the "optimal purifier" (optimal for question 
above) are defined in Section m and their fidelities are computed. The proof of 
the optimality claims is given in Section o. Finally, in Section 0, we determine 
the asymptotic behaviour for the optimal purifier, and the output rates. 

2 Figures of Merit and Main Results 

In this section we will state the optimization problems for purifiers mathemati- 
cally. A device (not necessarily a purification procedure) taking N qubit systems 
as input and producing M output qubits is described mathematically by a trace 
preserving, completely positive linear map ("cp-map") 

which takes input density matrices to output density matrices. Equivalently, 
we may work in the Heisenberg picture, using the dual T of T*, the unital (i.e. 



T(I) = I) cp-map 

T : Bin®^) -^ Bin®^), 

which is related to T, by tr(r(X)p) = tr(xr,(/9)). Here H = C^ is the one 
qubit Hilbert space, B{-) is the space of ah (bounded) operators on the cor- 
responding Hilbert space and i3,( • ) denotes the space of trace class operators. 
Since diniH = 2 < oo, the spaces S»('H) and B{'H) are just the 2 x 2-matrices, 
but it is nevertheless helpful to keep track of the distinction between spaces of 
observables and spaces of states. 

"Good purifiers" should make r»((i?,cr)'»^) very close to ct'^*=^. A simple 
figure of merit is the fidelity of the output with respect to the desired state in 
the worst case, i.e., 

^an(T) - inf tr (a®^-^T,((i?,a)®^))) , (2) 

where the infimum is over all one-particle pure states a. Similarly, we could 
pick any one of the outputs, say the one with number i, I < i < M , and test its 
fidelity. The worst case then gives the fidelity 

Tonc{T) = inf inf tr (a«r,((i?,a)®^))) , (3) 

where cr''-* = I(X>---®cr(X>---®I denotes the tensor products with (M — 1) 
factors "I" and one factor a at the i^^ position. We seek to maximize these 
numbers by judicious choice of T. Let us denote the optimal values by 

Tp''{N,M) = snp TiiT), (4) 

T 

where \^— "all" or \^~ "one" , and the supremum is over all unital cp-maps T with 
the specified number of inputs and outputs. 

For devices with variable numbers of outputs all these quantities become 
random variables, as well. Typically, one will seek to optimize the mean fidelity. 
It is then natural not to take the infimum in Equation (j^, but the mean. The 
case where no output is produced at all, is interpreted here as one output qubit 
in the completely mixed state. The resulting mean fidelity ||^ can be thought 
of as the fidelity ^ono(T) of a modified device T, which uses T, followed by a 
random selection of one of the outputs. Therefore, the problem of maximizing 
mean fidelity is exactly the same as maximizing J-onc{T) for devices with fixed 
output number M — 1, with optimal value jrP^'*^(iV, M). 

Rather than looking at the mean of the fidelity distribution of a device with 
variable number of outputs we could also look at its maximum. This corresponds 
to the problem in item n^ of the previous section. More precisely, one should omit 
the "worst case" infimum with respect to i in this case, and allow the device 
to either pick one of its outputs, or to declare failure. This leads to a device 
with only the two output numbers and 1, and the functional to be optimized 
is the fidelity of the "l"-output. We will denote the optimum for this problem 
by !Fa^^^^{N,0), with a slight abuse of notation expressing that this is the case 
with no demands on output numbers at all. 

It is clear that J-t'^^^{N, M) is a decreasing function of M, and that therefore 
the limit 

J^r^'^(iV, oo) = lim J^'°^'^(A^, M) (5) 



However, for tt=one, it is an interesting 

» CX). 



exists. For tt=all, this limit is zero 
quantity, which even goes to 1 as TV 

The results for the quantities T^Ti^, 0), T^T{N, 1), and T^'^iN, oo) are 
shown in Figure |l|. Of course, all these quantities also depend on the parameter 
describing the noise, which we have suppressed for notational convenience. It is 
fixed in the following graphs as A = 0.5 (resp. (3 — 0.549, see Section y|). It is 
clear that T^'^^^{N,M) -^ 1 for any N and M, as the noise level goes to zero 
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Figure 1: The three basic fidelities for the one-particle figure of merit: 
top: .Fr/(A^,0), middle: r^^{N,l), bottom: .Fr/(iV,oo) 



The leading asymptotic behaviour (as N — > oo) is of the form 

TZ^{N,M) CX l-^ + --- 

Co = (1-A)/A 

ci - (1-A)/A2 

Coo - {X+l)/X\ 



(6) 

(7) 
(8) 
(9) 



From these asymptotic results, a simple estimate for the all-particle fidelity 
criteria can be obtained: By Equation (^), 1 — J-'aii(T) < M{\ — J-onc{T)), 
where M is the number of outputs. Hence, for sufficiently small rate M/N one 
achieves good fidelity, even for the all-particle test criterion: 1 — J^y^{N, M) < 
M(l - J^^^{N,M)) < M(l - T^''{N,oo)) « ^Coo. Of course, the second 
estimate is rather crude, and a refined version will be given in Section o. The 
argument does show, however, that one may expect optimal all-particle fidelity 
to become a function of the output rate. This function will be computed in 



Section 3.3): for every fi > 0, we find the limit 



2A^ 



*(a*) 



lim 

M/N ^ti 



T^r{N,M) = 



2A2-^/i(l 
2A2 



A) 



if u< A 



ii H> X. 



(10) 



The function <I> is continuous and satisfies $(0) — 1 and <I>(oo) = 0, so at small 
rates purification is near perfect, but becomes arbitrarily bad at too high rates. 



In Figure g <& is plotted with the noise parameter A going in steps of 0.1 from 
to 1. The dotted Hne describes the performance of the natural purifier (see 
Section 4.1), which operates with rate /x = A. 
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Figure 2: Asymptotic fidelity <&(/i) for the all-particle figure of merit (px|). 
Curve parameter: A = .1, .2, ..., 1; dotted line: natural purifier 



3 Decomposition theory 

Many arguments in this paper are based on group theory, in particular the 
decomposition of tensor products of irreducible representations of SU(2). In 
this section we will summarize the relevant results which are needed throughout 
the paper. 



3.1 Reduction to fully symmetric case 

There are two reason why group theory is useful for us. First of all the depo- 
larizing channel R producing the noise is "covariant" which means that it does 
not prefer any particular polarization direction (basis in the underlying Hilbert 
space 7i = C^), and second we are looking at a "universal" purification problem, 
i.e. the purification devices T we are looking for should work well on an arbi- 
trary unknown input state a. Therefore, it is natural to look at those T which 
are covariant as well: T should work in exactly the same way on any input. 
Carrying this idea further it should also be impossible to single out any one of 
the input and output channels. Mathematically, these "natural conditions" are 
stated as follows: 

Definition 3.1. A unital, cp-map T : B{n®^'') -^ B{H'^'^) is called fully sym- 
metric if it is U(2) covariant, i.e. 

T{U®^AU*'^^) = u®^T{A)U*®^ VA G B{H®^'^) VC/ e U(2) 

and permutation invariant, i.e. 

T{t]Atj*) = T{A) \/ri e Sm VA e B(n®^') 



TTiA)T* = T{A) Vr e Sw VA e B{H®^). 

Here t] G Sm, t G Sat denote permutations of M respectively N elements and 
at the same time the corresponding unitaries on B{TC^^^) and B{'H®^), i.e. 
r]{ipi (X) ■ • ■ (g) i^Ai) = '0,,(i) ■ • ■ ® ^J,^{M) . 

We could have made this condition part of our definition of a purifier, and 
restricted the discussion to fully symmetric operations from the outset. How- 
ever, we have chosen to take the heuristic arguments at the beginning of this 
section more seriously: the kind of "universality" described there is already 
embodied in the figures of merit of Section ^ so it becomes a mathematical 
question whether optimal purifiers are indeed fully symmetric or else symmetry 
is broken, and a non-symmetric purifier can outperform all symmetric ones. 

We now argue that the optimal devices (with respect to J\)ne and J^sdi) may 
be indeed by assumed to be fully symmetric. To make this precise, note that 
•^aii(T') and J-onciT) are infima over expressions which are linear in T, and hence 
concave functionals. Therefore, averaging over many T's with the same figure 
of merit produces a T at least as good. Clearly, for all permutations 77 G Sm, 
r G Sat and t/ e U(2), the purifier T'{X) = TU®^T{riU*®^ XU^^ri*)U*'>^^T* 
has the same figure of merit as T . By averaging over these parameters (with 
respect to the appropriate Haar measures) we thus find a purifier, which is at 
least as good as T and, in addition, fully symmetric. Similar arguments apply 
for purifiers with variable numbers of outputs (although one has to be more 
careful in defining figures of merit). Therefore, we will restrict our discussion to 
fully symmetric purifiers from now on. 

3.2 Decomposition of tensor products 

The reduction to fully symmetric purifiers allows the application of techniques 
from group theory (especially representation theory of SU(2)) which simplifies 
our problems significantly. Consider in particular the A^— fold tensor product 

SU(2) 3U^ ^i/2{U)®^ = U®^ G B[n®^), 

of the spin-1/2, or the "defining" representation SU(2) 3 [/ h^ ''^i/2{U) = U E 
B(TL). It decomposes into a direct sum of irreducible subrepresentations 

7:i/2iU)®^ = U®^ = 7r,([/)®I (11) 

sei[N] 



ith 



and 



sei[N] 



^ } N even 
} TV odd 




[12) 



Here tTj denotes the spin-s irreducible representation of SU(2), Hs its 2s + 1- 
dimensional representation space, which we will identify in the following with 
the symmetric tensor-product Tif^^ , i.e. the 2s-qubits Bose subspace, and ICn.s 
denotes a multiplicity space, which carries an appropriate representation of the 
symmetric group Sat. 

3.3 Decomposition of states 

Consider now a general qubit density matrix p, which in its eigenbasis can be 
written as (/3 > 0) 

P^^^ = 2^^) ^"P ('^ y) ^ l^f^ ( l-^ ) 
- tanh(/3)|V')(V| + (1 - tanh(/3))il, ^ = f J 

The parametrization of p in terms of the "pseudo-temperature" (3 is chosen here, 
because it is, as we will see soon, very useful for calculations. The relation to 
the form of p = i?,cr initially given in Equation (Q) is obviously 

A = tanh(/3). 

The A^-fold tensor product p®^ can be expressed as p(/3)®^ = (2 cosh(/3))-^ exp(2/3L3) 
where 

B(W®^) 9 La = i (^3 ® I«(^-^) + ■ • • + 1^(^-1) ® a^) . (13) 

denotes the 3-component of angular momentum in the representation 7r®,2 ■ 
In other words, the density matrices are just analytic continuations of group 
unitaries, or "SU(2)-rotations by an imaginary angle 2i/3". This reduces the 
decomposition of p(/3)®^ to the decomposition (O) of the tensor product repre- 
sentation. Of course, analytically continued group elements are not normalized 
as density operators. Extracting appropriate normalization factors the decom- 
position becomes 

p(/3)®^= WM{s)p.sm^ ^ , 



with 



and 



sinh((2.s + 1)/?) 
sinh(/3){2cosh(/3))^ 



'^n{s) = „,„i^,^^,^'_,.^,^,.jv dim/CAf,^, (14) 



sinh(/3) „„„^r,flr('') 

smh{{2s + l)f3) 



p^(/^)- ..„.::^:...^ exp(2/34-^)■ 



(s) 

Here Lg denotes again the 3~component of angular momentum, now in the 
representation tt^. 

The PsiP) are normalized, i.e. trps(/3) = 1. Hence J^s'^Nis) = 1 and 
< WAr(s) < 1 due to the normalization of p(/3)®^. Together with the fact that 



the multiplicities dim ICn,s are independent of (3 we can extract from Equation 
(O) a generating functional for dimA^^jv^: 

2sinh(/3)(2cosh(/3))^ = 2 ^ sinh((2s + l)/3) dim/C^.s 

sei[N] 

= {e''-e-^){e^ + e-^f^ ^ (e^^^+i)'' - g-^^^+i)^) dim^^,„ 

sGl[N] 

obtaining 

dim/Cw s = 



N/2 + S + 1 \N/2 - sj ' 

provided N/2 — s is integer, and zero otherwise. The same result can be derived 
using representation theory of the symmetric group; see [Q, where the more 
general case dimTi = rf G N is studied. 

3.4 Decomposition of operations and optimal cloning 

Let us come back now to fully symmetric cp-maps T : B{Tt'^^^) ~* S(7i®^). 



Using the results of Subsection 3.2 it is easy to see that T can be decomposed 
into a direct sum 

T{A)^ Ts{A)®t (15) 

sei[N\ 

where the Ts : B{H®^) — >• B{'Hs) are unital cp-maps which are again fully 



symmetric (usi ng an obvious modification of Definition 3.1). Identifying, as in 
Subsection |3.2| , the representation space Tis with the 2s-fold symmetric tensor 
product Ti,^^ leads to the significantly simpler problem of decomposing fully 
symmetric, unital cp-maps Q : Biji,®^) -^ B{7i® ), which is already solved in 
[0. Hence we will state only the corresponding results here. In particular we 
have the following theorem: 



Theorem 3.1. Consider again the S-components of angular momentum L^ and 
Lg in the representations 7rr,2 respectively lis (cf. Subsection \3.S^ . 

1. For each fully symmetric cp-map Q : B{TC®'^^) -^ B{TLi_ ") there is a 
constant uj{Q) G M^ with Q{L^) ~ uj{Q)L^ . 

2. For each 2s G No there is exactly one fully symmetric Q2s with 

M 

— for 2s > M 

2s 
w((52s) = maxa;(Q) = ■( (16) 

-— ^ for 2s < M, 

2s + 2 ■' 

where the maximum is taken over the set of all fully symmetric cp-maps 
Q:B{n®^)^B{nf'). 



3. If M > 2s holds Q2s is given in terms of its pre-dual Qas* : B^,{Ti^ '*) — > 
B,{H®^) by 

Q2sM = §^^M(e ® I«(*^-2^))5m (17) 

where Sm is the projector from 7i®*^ onto the Bose subspace 7vf_ . 
4- For M < 2s the map Q2s is given by 

or in terms of its predual 

Q2s*i0)^tr2s-M0, (18) 

where tr2s-M denotes the partial trace over the first 2s ~M tensor factors. 

Note that the family of cp-maps Q2S defined in Equation (O) respectively 
(Iq) plays a very special role not only mathematically: Q2S describes the optimal 
way to increase (|l7| ) or dercrease (|lg) the number of qubits. More precisely Q2S* 
maps a finite number 2s of qubits in the same unknown pure state a to the best 
possible approximation Q2s*{o'^'^^) of the product state tr*^^^. The quality of 
Q2s*{<^^^'') is measured here by the fidelities 

aan(Q):=mftr(a®^^g,(a®2-)) 



or 



GoneiQ) :=mf inftr (a«g,(a®2^) 



If 2s > M holds (item Q) we simply have to discard 2s — M qubits to get exactly 
Q2s*(o''^^*) = (T®*^. If the number of qubits should be increased, i.e. M > 2s 
holds (item g), the target state cr*^*^ can not be reached. In this case Q2S is the 
optimal quantum cloning device described in ||l|, 3. 

4 Natural and optimal purifiers 

In this section we will introduce a particular class of purification maps which 
arise very naturally from the group theoretical discussion of the last section and 
which maximize, as we will see in Section o, the fidelities .7-aii and ^one- 

4.1 The definitions 

As a first step let us reinterpret the decomposition of p{l3)®^ discussed in Sub- 
section |3.3| in terms of the of cp-map 



B{'H%^')3 A, =: A ^ r""*(A) := Tf^\As) := 
se/[Af] sei[N] sei[N] 

:- As®le B{ns®lCN^s)=B{n''^''). (19) 

se/[JV] se/[JV] 
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Its predual maps the density matrix p[fi)®^ to 0sg/r^i wn{s)ps{P)- The lat- 
ter should be interpreted as a (normal) state on the von Neumann algebra 
®se/fJVl S{Hf'^'^). Hence T"^* is an instrument which produces with probabil- 
ity wn{s) the 2s-qubit state Ps{P) from the input state p{/3)^^ . This implies 
in particular that the number of output systems of T"'** is not a fixed parameter 
but an observable. We will sec soon that the fidelities of the output states Ps{P) 
are bigger than those of the input state p(/3)®^ provided s > holds. Hence 
we will call T"'** the natural purifier. 

The most obvious way to construct a device which produces always the same 
number of output systems is the composition of T"'** with the cloning operation 

B{n'^^)3A^Q{A)= Q2s(.A)e BiHf'). 

s£llN] sellN] 



Here the Q2s are the operations introduced in Theorem 3.1. Combining T"*^* 
with Q we get an operation 

Bin®"'') 3 A^ T°P'(yl) := (r"^*Q)(A) e B{n®^) (20) 

which produces, as stated, a fixed number M of output systems from N input 
qubits. Physically we can interpret T°^^{A) in the following way: First we apply 
the natural purifier to the input state p{l3)'®^ and we get 2s output systems 
in the common state Ps{P). If 2s > M we throw away M — 2s qubits and end 
up with a number of M. If 2s < M we have to invoke the 2s — * M optimal 
doner to reach the required number of M output systems. Although this cloning 
process is wasteful we will see soon that the fidelities JF^(T°p*) of the output 
state produced by T°p' are even the best fidelities we can get for any A^ — > M 
purifier. Hence we will call r°P* therefore the optimal purifier. 

4.2 The one qubit fidelity 

Now we will calculate the one qubit fidelity ^onc- Due to covariance of the 
depolarizing channel R the expressions under the infima defining Jo„(,(T) (and 
Ts\\{T)) in Equation \^ and (g|) depend for any fully symmetric purifier not on 
a and i. I.e. we get with R^a — p{l3): 



^an(T) = tr [a«^^T,(p(/3)®~)] and T,^,{T) - tr a(i)T,(p(/3) 



^®JV^ 



(21) 



with a = \il]){ip\. In the case of Tone the situation is further simplified by the 
introduction of the hlack cow parameter (cf. |l||) ^{6) which is defined for each 
density matrix on T^"**^ by 

7(e) = ^tr(2i30). 

To derive the relation of 7 to Tone note that full symmetry of T implies equiv- 
alently to (|l]) 



lillf^^MPf 



^_(r)=tr |-^V-.WlT^/..m®^^ 
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Since cr = (I + cr^)/2 holds with the Pauh matrix 0-3 we get together with the 
definition of Lj, in Equation (13) 

1 



''one I,-' ) 



1 



7[r*(p(/3)^^)] 



(22) 



In other words it is sufficient to calculate "(\T*{p{P)'^^)\ (which is simpler 
because SU(2) representation theory is more directly applicable) instead of 

'^onc v-^ ) • 

Another advantage of 7 is its close relation to the parameter A = tanh(/3) 
defining the operation _R, in Equation (||). In fact we have 

7(p(/3)^^) = 4tr(2i3p(/3)^^) = ^Niv{a^p{P)) = tanh(/3) = A. 



N 



N 



are 



In other words the one particle restrictions of the output state T(^p{f3)^^) 
given by 

7[r(p(/3)«^)]a+[l-7[np(/3)n]]^- 

This implies that j[T{p{f3)'^'^)] > A should hold if T is really a purifier. 

Let us consider now the natural purifier T"^' . Since the number of output 
qubits is not constant in this case we have to consider for each s € I[N^] the 
quantity J^oneiTs""^) (see Equation (|l^) for the definition of the T"'*') instead of 
one fixed parameter ^one(T'"'**) (in other words: The fidelity of T"^' is, as the 
number of output qubits, not a constant but an observable). According to the 
discussion above we get 



7(p.(/3))-^tr(24^V.(/3)) = 



1 tr(2L^''^exp(2/34"-')) 
2s tr(exp(2/3L^"^)) 



1 d 



1 d 



= :r:7;5l^tr(exp(2/34^^)) = -— (lnsinh((2s + l)/3) - lnsinh/3) 



2sdf3 

2s + 1 

2s 



coth((2s 



2s dp 
l)/5)-^coth/3 



(23) 



and hence 



*^ one V-^ s 



1 r 

2 

1 

2 



1 



1 



7(Ps(/3)®^ 
2s + 1 



2s 



coth((2s + 1)13) coth/3 



If s = 1/2 we have j{ps{P)) = tanh(/3) = A hence the (perturbed) input state 
p(/3) is reproduced. Taking the derivative with respect to s shows in addition 
that 7(ps(/3)) is strictly increasing in s. Hence T"^' really purifies (according 
to the remark above) and the best result we get if s is maximal. In the limit 
s ^ we find 7(ps(/3)) — which is reasonable because T"*** does not produce 
any output at all in this case {dimHs = 1 for s = 0). 

Let us apply these results to the optimal purifier. According to the definition 
of T°P' and T"^^* in Equations (|^ and {^ the decomposition of r°P* given in 
LSI) has the form 



T°P\A) = T'^\Q{A)) = Y. Q2.(A)®I= Y. T°P\A)®I, 



(24) 



se/[Af] 



se/[Ar] 
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hence T°p\A) = Q2siA). Together with @ we get 



^one(r°P*) = 



1+ J2 WNis)j[TTiPsW)] 

sei[N] 

1+ Y. WN{sh[Q2s*{Psim] 

seI[N] 

Y ■WN{s)fonc{M,(3,s), 

eilN] 



(25) 



where we have introduced the abbreviation 



fonc{M,P,s):=-[l+^[Q2siPsif3))] 



Together with Theorem |3.1| this imphes: 

2/onc(M, P,s)-1 = 7[Q2s*(P.(/3))] = ^ tr[2Q2siL3)ps{P)] 



^tr[2Li^^pM]- 



L0{Q2s)2s 

M 



i[psm- 



Inserting the values of uj{Q2s) and ^[ps{P)] from Equations dll) and ( p^ ) we 
get 



2/one(M,/3,s)-l 

■2s + l 



coth((2s + l)/3) - — coth^ 
2s 2s 



for 2s > M 



1 Af + 2 



('(2s + 1) coth((2s + 1)0) - coth/3] for 2s < M. 



(26) 



2s + 2 M 
Hence we have proved the following proposition. 

Proposition 4.1. The one-qubit fidelity J-^onciT"^^) of the optimal purifier is 
given by 



J^oneiT^P')^ Y, WN{s)fonciM,(3,s) 
si£l[N] 



(27) 



with /onc(-^, /?, s) from Equation (fe 



Note in particular that in the case M = 1 the one-qubit fidelity coincides 
with the expectation value of the fidelity of T"'^' in the state T^'^*(p(/3)®^) - the 
mean fidelity. Hence we can reinterpret the natural purifier as a device which 
produces exactly one output system (cf. ||^). 

4.3 The all qubit fidelity 

As in the one-qubit case the all-qubit fidelity of T'^^^ is an observable rather 
than a fixed parameter. Hence we have to calculate J^ediiT^'^*') for each fixed s. 
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Applying again Equation (21) we get 



^.n(T. )-tr(a p^i(3)) . —-^^-——^e 

g(2s+l)/3 _ g(2s-l)/3 I _ g-2;3 



g(2s+l)/3 _ g(2s+l)^ 1 _ g-(4s+2);3- 

Using the decomposition of r°P* given in Equation (EJ) we get for the opti- 
mal purifier something similar as in the last subsection: 

^an(T°P')- ^ WN{s)tT[a^'^'T°t{ps{m (28) 

sei[N] 

= ^ WNis)tT[a^''Q2s.{Psm) 
sei[N] 

However the calculation of 

/an(M,/3,s):=tr[(T«*^Q2s*(ps(/3)) 

is now more difficult, since the knowledge of Q2s{L3) — uj{Q2s)L^ is not suffi- 
cient in this case. Hence we have to use the explicit form of Q2s in Equation 
(0) and (|l|). For 2s < M this leads to 



/aii(M,/3,s) 



{i:^'',SM{Ps^l^^''-''hsMr'') 



2s + l 
A/ + 1 
2s -1-1 
M + 1 
2s + 1 1 - e-2/3 



(V'®'^,(p,®I^(^^-2^))V'®''') 



2s -Hi 
M + 1 



{^p^'^^Ps^n 



M + 1 l-e-(4^+2)/3- 
For M < 2s we have to calculate 

/an(s,M,/3) = tr [a®*^g2s*(ps(/3)) 



= tr 



i2s 



K"0p.(/3) 



= tr 



p,(/?)(5m[(|V'®^'')(V'^^''|)®: 



i{2s-M)- 



M 



(29) 



We will compute the operator Q2s{o'^'^^) in occupation number representation. 
By definition, the basis vector "|n}" of the occupation number basis is the 
normalized version of Sm"^, where ^ is a tensor product of n factors "0 a-nd 
(M — n) factors 0, where 4> = d) denotes obviously the second basis vector. 
The normalization factor is easily computed to be 



M 



-1/2 



5m(V'''" ® 0»(^^-")) = ( '-^^ ) \n). (30) 



We can now expand the "I" in Equation (g9|) in product basis, and apply (|30|), 
to find 



SMm'^^'){c^®^'\) ® i®("^-'')]5m = Y. 



K 



2s - M\ /2s 



K -M \K 



\K){K\. 
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Now ia is diagonal in this basis, with eigenvalues rriK = [K — s), K = 
0, . . . , (2s). With p,{P) from (P) we get 



/an(M,/3,s) 



Together with 



-2/3 



■E 



1 _ e-(4s+2)/3 A^ \K -M l\K 

K 



2s - M\ /2s 



^2f3{K-s) 



for M < 2s. 



2s - M\ (2s 



^K -M) \K 
we get 



(2s -M)! K\{2s-K)\ 

[K - M)\{2s - K)\ (2^)! 



/aii(M,/3,s) 



1 - e-^'^ /2s 
1 - e-(4«+2)/5 U/ 



E 



e 



^2/3(/f-s) 



Summarizing these calculations we get the following proposition: 

Proposition 4.2. The all-quhit fidelity J-a.\\{T°^^) of the optimal purifier is 
given by 



.Fall(T°P*)= Y. WN{s)fonciM,P,s) 

sei[N] 



(31) 



where /au(M, /3, s) is given by 



2s + 1 1 



-2/3 



/an(Af,/3,s) = <^ 



M + 1 1 



-(4s+2)/3 



1 



= -2/3 



2s 



1 _ g-(4.+2)/3 U/ 



E 

K 



M < 2s 



^2PiK-s) M>2s. 



(32) 



5 Solution of the optimization problems 

Now we are going to prove the following theorem: 

Theorem 5.1. The purifier T°^^ maximizes the fidelities J-onc{T) and !Fa\\{T). 
Hence the optimal fidelities T^^^{N,M) and T^^''{N,M) defined in Section^ 
are given by Equation l[2\) and (|j'j|j. 



Proof. Note first that the funtionals J-"onc and TaW are, as infima over continuous 
functions, upper semicontinuous. Together with the compactness of the set of 
admissible T this implies that the suprema J-4^^^{N, M) from Equation (H) are 
attained. In other words: optimal purifier T with J-^{T) = !F^^'^{N,M) exist, 
and we can assume without loss of generality that they arc fully symmetric 
(according to the discussion in Section 3.1). Hence we can apply Equation (Uw 
and the decomposition (p^ to get in analogy to (|2^) and (p8[) 



'^ onct,-^ J 



1+ Y. WN(sh[Ts*{psm] 

seilN] 
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and 



^an(T)= V WN{s)tT[a'^^'Ts4psil3))]. (33) 



sei[N] 



The last two Equations show that we have to optimize each component Ts 
of the purifier T independently. In the one qubit case this is very easy, be- 
cause we can use Theorem 3J_ to get TsiL^) = io{Ts)L]^^ and 7 [Ts*(/9s (/?))] = 



uj{Ts)tr(^L]^ Ps{P)). Hence maximizing 7[Ts»(ps(/?))] is equivalent to maximiz- 
ing uj{Ts). But we have according to Theorem 



msLXUi{Ts) = uj{Q2s) = 



T 




2(s + l) 



for 2s > M 
for 2s < M, 



which shows that T^TiN, M) = J^onc(r°P') holds as stated. 

For the many qubit-test version the proof is slightly more difficult. However 
as in the ^one-case we can solve the optimization problem for each summand in 
Equation (p3) separately. First of all this means that we can assume without 
loss of generality that Ts* takes its values in B{7if ) because the functional 

/,(T,):=tr(a®^'^r,,(p,(/3))) (34) 

which we have to maximize, depends only on this part of the operation. Full 
symmetry implies in addition that Ts,,{ps{P)) is diagonal in occupation number 
basis (see Equation (|3^)), because TsripsiP)) commutes with each ns'{U) (s' = 
M/2, U e U(2)) if TT,{U) commutes with p^(/3). 

If M > 2s this means we have T'sH,(yOs(/3)) = K»cr®*^ + r, where r, is a 
positive operator with cr®^^r, = r^.a®'^ = 0. Inserting this into (|3J) we see that 
faiTg) = K*. Hence we have to maximize k,. The first step is an upper bound 
which we get from the fact that tr(cr®*^ps(/3))l — Ps{P) is a positive operator. 
Since rs,(I) = (2s -I- 1)/{M + 1)1 (another consequence of full symmetry) we 
have 



< r(tr(a®2v^(^))l _ p^(/3)) = |i±i- tr(a«^V.(/3))l - ^<y^'' 



Multiplying this Equation with (j®*^ and taking the trace we get 

^^*<§-^tr(a®^V(/?)). (35) 

However calculating fs{T°P^) we see that this upper bound is achieved, in other 
words T°P* maximizes fs- 

If M < 2s holds we have to use slightly different arguments because the 
estimate ( pq ) is to weak in this case. However we can consider in Equation (B3) 
the dual T, instead of Ts* and use then similar arguments. In fact for each 
covariant Ts the quantity ^5(0-®*'^) is, due to the same reasons as Ts^,{ps{(3)) 
diagonal in the occupation number basis and we get Ts{(J®^'^) = na'®'^^ + r 
where r is again a positive operator with r = X^ntJo ^ri\n) (|n) denotes again 
the occupation number basis) and k is a positive constant. Since Ts is unital 
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we get from I — a'^^ > the estimate < k < 1 in the same way as Equation 
(Bq). Calculating T°^^{a®-^) shows again that the upper bound k = 1 is indeed 
achieved, however it is now not clear whether maximizing k is equivalent to 
maximizing fs(Ts). 

Hence let us show first that k = 1 is necessary for fs{Ts) to be maximal. 
This follows basically from the fact that T, is, up to a multiplicative constant, 
trace preserving. In fact we have 

tr(r,(a®*'0) - tr(T,(cT«*Ol) - tr(a«*^T„(I)) = ^J±i. 

This means especially that k + tr(r) = (2s + l)/(Af + 1) holds, i.e. decreasing 
K by < e < 1 is equivalent to increasing tr(r) by the same e. Taking into 
account that Ps{P) = X]n=o ^nl") holds with /i„ = exp(2/3(n — s)), we see that 
reducing k by e reduces fs{Ts) at least by 

e{tr{<j'^'^pM) - tr(|2s - l)p.(/3))) = e{e'^^ - e'^-D/^) > 0. 

Therefore k = 1 is necessary. 

The last question we have to answer, is how the rest term r has to be chosen, 
for /s(T's) to be maximal. To this end let us consider the slightly modified fidelity 
f^{Ts) = tr(7Ucr®^)cr'»2s-) (which is in fact related to optimal cloning; see 
and Section |3.4| ). It is in contrast to /s(Ts) maximized iff n — 1. However the 
operation which maximizes f{Ts) is obviously the optimal M —> 2s doner (up 
to normalization) which is according to unique. This implies that k = 1 fixes 
Ts already. Together with the facts that k = 1 is necessary for fs{Ts) to be 
maximal and k = 1 is realized for T°p* we conclude that raax fs{Ts) = /s(r°P*) 
holds, which proves the assertion. D 



6 Asymptotic behaviour 

Now we want to analyze the rate with which nearly perfect purified qubits can be 
produced in the limit N —> oo. To this end we have to compute the asymptotic 
behaviour of various expectations involving s. It turns out that it is much better 
not to do work with the explicit expressions of these expectations, as sums over 
expressions with many binomial coefficients, but to go back to the definition, 
and use general properties of expectations of p^^ . This has the added advantage 
of being easily generalized to Hilbert space dimensions d > 2, so we expect the 
method to be useful in its own right. We collect the basic statements in the 
following subsection, applying them to the concrete expressions in subsequent 



6.1 Convergence of weights to a point measure 

In the classical case the general theory alluded to above is nothing but the the- 
ory of asymptotic distributions for independent identically distributed random 
variables (Laws of large numbers of various sorts). In the quantum case this 
theory has been developed in the context of the statistical mechanics of general 
mean- field systems H . Of this theory we need only the simplest aspects (con- 
vergence to a point measure), and not the more advanced "Large Deviation" 



17 



parts, in which it is shown how the probabiHty of deviations from the hniit 
decrease exponentially fast. 

Consider operators of the form A^ = (1/^) X]i=i '^^*^ where a*^'^ denotes 
the copies of a fixed operator on Ti, acting in the «*^ tensor factor of Ti,^^ . 
It is clear that the expectations ti{p'^^ An) = tr(pa) are independent of N. 
Now consider products of a finite number of such operators and expand the 
expectation into the average over all terms of the form iT:[p®^ a^^'>b'^^^ c'^^^ " ■ ■ )• 
It is easy to see that for large N the majority of these terms will be such that all 
indices i, j, /c, . . . are different, and for such terms the above expression is equal 
to tr(pa) tr(/96) tr(pc) ■ • ■ . So this will be the limit of the expectation of the 
product A^BnCn ■ ■ ■ as N ^ co (for precise combinatorial estimates, see Q). 
Of course, this allows us to compute the asymptotic expectations for arbitrary 
polynomials, and by taking suitable limits of arbitrary continuous functions of 
Hermitian operators. There is an abstract non-commutative functional calculus 
describing exactly these possibilities (see appendix of g). However, for our 
purposes it is sufficient to say that all combinations of algebraic operations and 
continuous functions of a Hermitian variable (evaluated in the usual spectral 
functional calculus) are in this class. 

For the case at hand, note that the angular momentum operators Lk as in 
Equation ( |l3| ) are of the form NAjq therefore, for any sequence of functions /at 
of three non-commuting arguments (this means that in writing out /tv we have 
to keep track of operator ordering), which converges to a limit function, /oo, we 
get 

Jmi^tr(p«^/^(^,^,f))=/^(tr(p^),tr(p^),tr(p^)). (36) 

Note that the function foo is just evaluated on numbers (operators on a one- 
dimensional space) so all operator ordering problems disappear in the limit. 
This is the huge simplification which makes mean-field theory so accessible. 
The limit formula will be applied to functions of "2s" , the number of outputs 
from the natural purifier, which can itself be written as a function of this sort. 
It is, of course, constant on each summand of the decomposition (O), so it is a 
function of the Casimir operator L" = s{s + 1): 



'^-^=9N{^,^,^)^^^{L/Nf+N-^-l/N 



ffoo(a;i,a;2,a:3)- lim ^^{xf + N-^ - l/N = 2\x\ 



N 



9oo (tr(p^), tr(p^), tr(p^)) = g^iO, 0, A/2) = A = tanh/3, (37) 

when p — p{(3) is given by eq.([l^). Functions of g then also lie in the relevant 
functional calculus, so we get the following statement, taylored to our need in 
the following subsections. In it we have already encorporated further, straight- 
forward approximation arguments, using uniformly convergent sequences of con- 
tinuous functions to establish upper and lower bounds separately. 

Lemma 6.1. Let /n : (0, 1) — > M, A^ G N 6e a uniformly bounded sequence of 
continuous functions, converging uniformly on a neighborhood of X = tr(/9(/3)CT3) 
to a continuous function f^o, and let wn{s) denote the weights in Equation ^ 
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Then 



Km V zi;jv(s)/jv(2s/7V) = /00(A). 

N—^OO ^ — ^ 

sei[N] 



(38) 



In the language of measure theory this is saying that the probabihty measures 
X^s fi'Ar(s)(5(a; — 2s/N)dx on the interval [0,1] converge to the point measure 
5{x — X)dx. Graphically, this is shown in Figure 
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Figure 3: Convergence of wn{s) to a point measure (A = .b, N — 10, 100, 1000). 
Discrete points joined, and rescaled for total area 1 



6.2 The one particle test 

Let us analyze first the behaviour of the optimal one-qubit fidelity J-^^{N, M) 
in the limit M -^ cxo. Obviously only the M > 2s case of fonc{M, /?, s) is relevant 
in this situation and we get, together with Equation (CT), the expression 



sei[N] 



1 



2s + 2 



(2s + l)coth((2s + l)/3) -coth/3 



which obviously takes its values between and 1. To take the limit A^ ^ 00 we 
can write 



N 



hm ^-r(7V,oo)= lim Y. y^N{s)fN,o.ijj) 

s6/[Af] 



with 



/jv,oo(a;) = - 



1 



Nx^2 



( A^a; + 1 ) coth ( (iVx + 1 )/?) - coth /3 



The functions Jn,oo are continuous, bounded and converge on each interval (e , 1) 
with < e < 1 uniformly to /oo.oo = 1- Hence the assumptions of Lemma 3.1 
are fulfilled and we get 

hm .f--r(7V,oo)==/oo,oo(A) = l 
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as already stated in Section 0. This means that we can produce arbitrarily good 
purified qubits at infinite rate if we have enough input systems. 

To analyze how fast the quantity ^^^{N, cx)) approaches 1 as A^ ^ oo let 
us consider the limit 



2s, 



hm 7V(l-.F-r(7V,^))= E WNis)fj,M-)^^ 



Af— >oo 



seilN] 



(39) 



with /jv.oo = -^(1 — /at.oo)- The existence of this limit is equivalent to the 
asymptotic formula 

where, as usual, o (-i-) stands for terms going to zero faster than -^. Lemma |6.l| 
leads to Coo/2 = foc,ooW with foo,oo = liniAr^oo /at^oo uniformly on (e, 1). To 
calculate /oo,oo note that 



In, ooix) 



N 



Nx- 



N coth 13 

Nx + 2 



Rest 



holds, where "Rest" is a term which vanishes exponentially fast as A^ -^ oo. 
Hence with coth/3 — 1/A we get 



2/oo,oo(A) 



1 + A 
A2 



The asymptotic behaviour of !F^'^{N, 1) can be analyzed in the same way. 
The only difference is that we have to consider now the 1 = M < 2s branch of 
Equation (Eq). In analogy to Equation (p9) we have to look at 



,2s. 



Cl 



hm 7V(1 - .^°r(7V, 1)) = J2 ^Nis)fNA-)= „ 



sei[N] 



with fpf,! = N{1 - fN.i) and 



fNA{x) = - 



1 \(Nx + 1) coth((A^a; + l)f3) - coth/3l 

Nx '- \ / J 



For /oo,i we get 



^-^(^) = 2(v + ;a)- 



Using again Lemma 3.1 leads to 



(40) 



Cl = 2/oo,i(A) = . 

Finally let us consider J^^^^{N,0). Here the situation is easier than in the 
other cases because ^™|'^(iV, 0) equals the fidelity of the best possible output 
of the natural purifier, i.e. 



•^nr(A^,0) = i 



1 - ^ [(A^ + 1) coth((A^ + l)/3) - coth/3] 



/^.i(l)- 
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Hence we only need the asymptotic behaviour of /jv,i(a^) at x = 1. Using 
Equation (|^ we get 



^r/(7v,o) = i 



l-A 1 



A 2N 
This concludes the proof of Equations (^ to (Q). 

6.3 The many particle test 

Consider now the many-qubit fidelity J-aii- Although, like J^onc, it lies between 
zero and one, and would attain the value 1 precisely for a (non-existent) ideal 
purifier, both quantities behave quite differently, when we use them to compare 
states in systems of varying size. We are looking here at the two kinds of 
fidelities for an M-particle output state pm with respect to a one-particle pure 
state given by the vector -0, namely 

Fan = (^®^^PM^®*^)-trpM(|^^*0(^®*'l) ,and 
F, = (7/;,pi^V)=trpM(l®---(|V'>(^l)^®l) , 

where p\^ denotes the restriction of pu to the ?'^ tensor factor. Let paii and 
Pi denote the projections whose pM-expectations appear on the right hand side 
of these Equations. These projections commute, and paii is the intersection (in 
the commuting case: the product) of the pi in the lattice of projections. This 
corresponds to the union of the respective complements, i.e., 

i 

Taking expectations with respect to pM, we find that supj (1 — Fi) < (1 — Fail) < 
^,■(1 — Fi) < Msupj(l — Fi). For the two figures of merit introduced in Section^ 
this implies 

(1 - ToneiT)) < (1 - T^niT)) < M(l - .Fonc(F)) , (41) 

for every purifying device T. Hence, for fixed N the two figures of merit are 
equivalent to within a factor . But the upper bound becomes meaningless in the 
limit M — > 00, so it is not clear at all whether we can bring the fidelity JF^ji(T) 
close to one for an increasing number of outputs. 

As a consequence of this analysis it is necessary to perform the limit TV, M — > 
oo more carefully as in the one qubit case. We will consider therefore the limits 
N -^ oo and M — > oo simultaneously, while the quotient M/N approaches a 
constant p,, i.e. we will calculate the function <&(/i) defined in Equation (M). 
The first step in this context is the following lemma, which allows us to handle 
the {l]y^ J2k (m)^'"^''^'"^ term in Equation ^). 
Lemma 6.2. For integers M < K and z G C, define 

^ ^ R=M ^ ^ 
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Then, for \z\ < I, and c > 1; 

lim ^{K, M, z) 



M,K^^ • ■ ■ ■ 1 - (1 - C)Z 



Proof. We substitute R^-^ {K — R) in the sum, and get 
$(i^, M,z) = Y^ c{K, M, R) 



z^. 



R=0 



where coefficients with M + R > K are defined to be zero. We can write the 
non-zero coefficients as 

c[K,M,R) . y [ m)= KMK R-My. 

_{K -M){K -M ^l) {K -M -R+l) 
~ K (K-1) (K-R+l) 

Since < c{K, M,R) < 1, for aU K, M, R, the series for different values of M, K 
are all dominated by the geometric series, and we can go to the limit termwise, 
for every R separately. In this limit we have M/{K — S) —> c for every S, and 
hence c{K, M, i?) ^ (1 — c)^. The limit series is again geometric, with quotient 
(1 — c)z and we get the result. D 

To calculate now $(//) recall that the weights wn{s) approach a point mea- 
sure in 2s/N =: x concentrated at A = tr{p{(3)a3). This means that in Equation 
( pl| ) only the term with 2s = XN survives the l imit . Hence if /i > A we get 



M > XN = 2s. Using Equation (|32D and Lemma 6.1 we get in this case 



We see that $(^) ^ for ^ ^ oo and ^{fj,) ^ 1 — exp(— 2/3) for /x ^- A. 

If < /i < A we get M < XN = 2s, w hich means we have to choose Equation 



2|) for /an(M, /3, s). With Lemma 6^ and Lemma |6.l| we get 

1 - e-2/3 



^^i) = 



1 - (1 - /i/A)e-2/3 



which approaches 1 if /i ^ and 1 — exp(— 2/?) if /i ^ A. Writing this in terms 
of A = tanh/3, we obtain Equation (|o|). 

6.4 Estimating the many particle fidelity in terms of one 
particle 

In Section g we motivated the observation that the the best all-particle fidelity 
is a function of the rate (and not identically equal to 1) by estimating the all- 
particle fidelity in terms of the one-particle fidelity. Since the latter quantity 
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tends to be more easily computable it is of some interest for further investiga- 
tions, how good that estimate actually is. The estimate mentioned in the text 
before Equation (^ amounts to 



$(m) > 1 - f Co 



1 



2A ' 



(42) 



However, the same basic estimate via Equation (Bl|) gives even more informa- 
tion: 



a>(M)>l- Jim M{1-TZT(N,M)) 



> 



1-fi lim V wn{s) N{l-fonc{f^N,l3,s)) 



se/[Ar] 



A^(l-A) 

1 ttt: II II < A 



.2A2 , 
/i(l + A) . 



(43) 



2A2 



if /i > A, 



where the evaluation of the limit was carried out with the same technique based 
on Lemma 6.1 used in the previous sections. Figure displays the lower bounds 
(H) and (||) together with the exact result (|l|). 



0.8 



0.6 - 




0.4 



0.2 - 



Eq. (10) 



Figure 4: The lower bounds (42) and ( [43|) together with the exact result ( |l(]| ) 
for the all-particle test fidelity as a function of the rate (A — .5) 



It is apparent that these bounds are rather weak, and in fact completely 
trivial for large rates. Hence all-particle fidelities contain new and independent 
information about purification processes, which is not already contained in their 
one-particle counterparts. 
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